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We obtain the non-collinear ground states of a triangular-lattice antiferromagnet with exchange 
interactions up to third nearest neighbors as a function of the single-ion anisotropy D. At a critical 
value of D, the collinear f f J.J, phase transforms into a complex non-collinear phase with odd-order 
harmonics of the fundamental ordering wavevector Q. The observed elastic peaks at 2-7TX — Q in 
both Al- and Ga- doped CuFeCh are explained by a "scalene" distortion of the triangular lattice 
produced by the repulsion of neighboring oxygen atoms. 

PACS numbers: 75.30.Ds, 75.50.Ee, 61.05.fg 



The non-collincar and multifcrroic ground states of 
frustrated magnetic systems continue to attract intense 
interest. Due to the strong coupling between the electric 
polarization and the non-collinear spin states, improper 
ferroelectric materials hold great technological promise 
[l|]. However, more than one physical mechanism may be 
responsible for their ferroelectric behavior. An electric 
polarization P perpendicular to both the spin rotation 
axis Si x Sj and the wavevector Q is predicted for ferro- 
electrics with easy-plane anisotropy and spiral spin states 
like i?Mn0 3 (R = Tb or Y) Q. But a simple spiral state 
is not possible for ferroelectrics based on materials with 
easy-axis anisotropy, like CuFeC>2 d H IH and MnW0 4 
;6j. For Al- or Ga- doped CuFeC>2, a modulation of the 
metal-ligand hybridization with the spin-orbit coupling 
7, 8 1 may produce the observed electric polarization P 
3 1 4 , [f| parallel to both the spin rotation axis and the 
wavevector. 

In order to clarify the nature of the ferroelectric cou- 
pling, it is essential to understand how the non-collincar 
ground state of an easy-axis ferroelectric evolves with 
doping. In this paper, we show that the non-collincar 
ground state of CuFeC>2 contains significant odd-order 
harmonics of the fundamental ordering wavevector Q w 
0.867TX Q. The observed elastic peaks at both Q and 
27rx Q f=a 1.14-7TX 0, [i| are explained by a distortion 
of the triangular lattice associated with the repulsion of 
neighboring oxygen atoms. 

Due to geometric frustration, simple antiferromagnetic 
(AF) order is not possible on a two-dimensional triangu- 
lar lattice with AF interactions J\ < between neigh- 
boring sites. When the easy-axis anisotropy D along the 
z axis is sufficiently large, however, the anisotropy en- 
ergy —DJ^iSiz 2 favors one of several collinear states. 
For classical spins, Takagi and Mekata [l(| demon- 
strated that the ft 1 1 state sketched in Fig.l is stable 
over the range of J2/IJ1I and J-s/\J\\ plotted in the in- 
set to Fig. 2 (a), where J 2 and J3 are the second- and 
third-neighbor interactions indicated in Fig.l and longer- 
ranged interactions are neglected. The fffj, phase with 
wavevector Qo = 7rx appears in pure CuFeCh 11, l3| for 



FIG. 1: (Color online) The exchange constants J{ , J2, and 

J3, and the oxygen displacements responsible for the scalene 

(3) 

distortion of the lattice. The thick bonds J{ — J\+Ki form 
a zig-zag pattern. Also shown is the ffff state with f (solid) 
and I (open) spins. For K\ > 0, j[ 3 ^ couples the same spin 
along each zig-zag. 



magnetic fields below about 7 T. 

With increasing Al concentration, the spin-waves 
(SWs) of CuFci_ K Al a; 02 soften on either side of Q at 
wavevectors Q± ~ Qo ± 0.14ttx [U, El- A similar SW 
softening occurs on a triangular lattice AF when D is 
lowered while the exchange constants are fixed [H} . For 
an Al concentration x greater than x c « 0.016 or an 
anisotropy D lower than D c w 0.3 1 J\ |, the ffff phase be- 
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comes unstable and a non-collinear phase appears [1 
with the dominant wavevector Q = Q_ f=a 0.867TX. 

Because the AF interactions between adjacent hexag- 
onal layers of CuFeC>2 are not frustrated, the essential 
physics of CuFe02 is captured by a two-dimensional 
triangular-lattice AF with energy 



E — — - JijSi - 



(1) 



A classical approximation for the 5 = 5/2 spins of the 
Fe 3+ ions incurs only small errors, so each spin Si = 
S(Ri) is treated classically with |S(R;)| = 1. We include 
exchange couplings Jy up to third nearest neighbors. 
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Monte-Carlo simulations were recently used to s tudy 
the complex non- collincar (CNC) phase of Eq.© [l8|. 
Those simulations indicated that the CNC phase inter- 
ceded between the ftll phase at high D and a spiral 
phase at small D. Several "trial" spin functions were 
constructed to minimize the energy of Eq.([T]), including 
functions with all three spin components. Those trial 
functions were motivated by the Fourier peaks in the 
Monte-Carlo solution [HI at wavevectors Q = 0.877TX, 
1.137rx = 27rx — Q, and 1.387rx f=a — 3Q + G where 
G = 4-7TX is a reciprocal lattice vector. 

The trial function with the lowest energy contains odd- 
order harmonics in an expansion of the spin: 

S Z (R) = a{cos(Qx) + J2 c 21+icos(Q(21 + 1)x) 

1=1 

+ J2 B 2i+icos(Q'(2l + l)x + ^)Y (2) 



i=o 



(3) 



where Q' = 2?rx - Q. Notice that S(R) = S(x) de- 
pends only on x. The anharmonic coefficients C2i+i>i 
reflect the deviation from a pure cycloid with S(x) = 
(0, sin(Qx), cos(Qa;)) ; the coefficients B 2 i+i are produced 
by a lattice distortion with period 1, as discussed further 
below. The amplitude A is fixed by the constraint that 
max|5' z (a;)| = 1 and the lattice constant is set to 1. Keep 
in mind that Q and Q' are distinct wavevectors not re- 
lated by a reciprocal lattice vector. 

Like the observed multiferroic phase the CNC 
phase of Eqs.((2|) and ([3]) is coplanar with the spin ro- 
tation axis S(x) x S(x + 1/2) parallel to Q along the x 
axis. A uniform rotation of S(x) about the z axis would 
not cost any anisotropy energy but would cost magneto- 
clastic energy due to the distortions discussed below. 

As shown in Ref. 18(, the dominant wavevector of the 
CNC phase coincides with the wavevector of the dom- 
inant SW instability of the phase. Depending 
on whether the exchange parameters { J%/\Ji\, Js/\Ji\} 
fall within regions 41 or 411 plotted in the inset to 
Fig. 2(a), the dominant SW instabilities occur at a vari- 
able wavevector Q (region 41) or at (47r/3)x (region 
411) The exchange parameters used in Fig. 2(a) 
(J 2 /\Ji\ = -0.20, J 3 /\Ji\ = -0.26) fall within region 
411; the exchange parameters used in Fig. 2(b) (J2/I Ji\ = 
-0.44, J 3 /\Ji\ = -0.57) fall within region 41. The lat- 
ter are believed to correspond approximately to the ra- 
tio of exchange parameters in pure CuFe02 [3]- With 
those parameters, the dominant SW instability of the 
tti-l phase and the dominant ordering wavevector of the 
CNC phase are both Q « Q 0.14ttx = 0.86?rx. The 
wavevector of the third harmonic 3Q is then equivalent 
to 1.42ttx 0. 

The classical energy E was minimized within a unit cell 
of length 5,000 with open boundary conditions in the x 



direction. Doubling the unit cell has no noticeable effect 
on the amplitudes plotted in Figs. 2 (a) and (b). In the 
absence of a lattice distortion, the amplitudes B21+1 arc 
negligible but the higher harmonics C2i+i>i are signif- 
icant. For all D/\Ji\, the trial spin configuration has a 
lower energy than the Monte-Carlo state. Notice that the 
anharmonicity in region 411 is much weaker than in re- 
gion 41. For the parameters of Fig. 2(a), only the third and 
fifth harmonics C3 and C5 are significant; for the param- 
eters of Fig.2(b), harmonics above C7 can be neglected. 
The S z component of the anharmonic CNC phase within 
region 41 is sketched in the inset to Fig. 2(b). This phase 
retains some of the Ising character of the till phase with 



0.72 at D = 0.3. 



Within region 411, Q = 47rx/3 depends on neither the 
exchange parameters nor the anisotropy; within region 
41, Q is relatively insensitive to D but sensitively de- 
pends on the ratio of exchange parameters, as discussed 
in Refs.[l8| and [19( 1 . For the parameters used in Fig. 2(b), 
Q w 0.857ttx. 

With decreasing D, the anharmonicity decreases in 
both regions 41 and 411. Pure spirals with €21+1^1 = 
and (S z 2 ) = 0.5 are recovered as D — > 0. So the phase 



diagram provided by Fig. 3 of Ref.[l8[ should be revised 
to eliminate the sharp boundary between the CNC phase 
and the spiral region. Comparing the numerical results 
obtained for the trial spin configuration of Eqs.© and 
Q with the earlier Monte-Carlo results [3] using the 
parameters of Fig. 2(b), we find that the critical value D c 
below which the ]][[ phase disappears increases from 
0.295|Ji| to 0.317| Ji|. When D/\.h\ = 0.1, the energy 
of the Monte-Carlo phase is E/N = — 1.284| Jx| whereas 
the energy of the anharmonic CNC phase is — 1.295 1 Jx|. 
These energies can be compared with the energy of the 
ttU phase in Fig.3(b). 

The observed 2-7TX — Q peak in the elastic neutron- 
scattering measurements 0, [f| with amplitude \Bi\ 2 is 
absent for a non-distorted lattice. This clastic peak re- 
quires a lattice distortion with a wavevector of q = 27rx 
or a period of 1. The most likely source of that distor- 
tion is the repulsion of neighboring oxygen atoms shown 
in Fig.l. For each pair of oxygen atoms, one lies below 
the hexagonal layer and the other lies above. The dis- 
placement of oxygen atoms pictured in Fig.l produces 
a "scalene" distortion of the triangular lattice that has 
been observed in both the low- field ttll phase of pure 
CuFeC>2 as well as in the field-induced multiferroic phase 
above 7 T (2Cj. It has also been reported in the ferro- 
electric phase of Al-doped CuFe02 [4[. The displacement 
expands the lattice in the x direction 2l| , which was ob- 



served in pure CuFe0 2 Hill. 



To study the distorted phase, we take J^' = = 

(3) 

J 1 — Ki/2 and Jf = Jx + K\, which ensures that the 
average bond strength remains J±. The energy K\ mea- 
sures the degree of distortion of the nearest-neighbor ex- 
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FIG. 2: (Color online) The harmonic amplitudes C21+1 as 
a function of D/\,h\ for (a) J 2 /\Ji\ = -0.20 and J 3 /\Ji\ = 
-0.26 or (b) J 2 /\Ji\ = -0.44 and J 3 /\Ji\ = -0.57. Inset 
in (a) is the phase diagram indicating regions 41 and 411 dis- 
cussed in the text with the red dot corresponding to the (a) 
parameters and the blue dot to the (b) parameters. In (b) 
we sketch the TTi-l phase stable above D/\J\ \ = 0.32 and the 
CNC phase [S z component only) at Dj\J\ \ = 0.3. 
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FIG. 3: (Color online) (a) The amplitudes B21+1 and C21+1 
versus the distortion K\ for D/\J\ \ = 0.1. Inset are the spin 
populations P(S Z ) for the spin states at Ki = (long dash) 
and Ki/\Ji\ — 0.05 (short dash) as well as their difference 
(solid), (b) The energy E/N\Jx\ for the predicted CNC phase 
(solid), the tTll phase (short dash) and the "proper helix" 
[3| (long dash) . The dot along the E axis is the energy of the 
Monte-Carlo simulation [l8l ]. Other parameters as in Fig. 2(b). 



change. When K\ > 0, the oxygen displacements weaken 
(3) 

the AF coupling J\ and strengthen the AF couplings 
J^ 1 ' and J-p . The bonds form a zig-zag pattern 
with wavevector 2-7TX or a period of 1. The lattice distor- 
tion lowers the energy of the Tti-L phase, which is given 
by E/N = Ji - J 2 + J 3 - D - 2K 1 for K x > and 
E/N = J 1 -J 2 + J 3 -D + Ki for K x < 0. For K x > 0, 
the | or I spins prefer to lie on the same zig-zag cou- 

(3) 

pled by the weakest AF bond J\_ , as shown in Fig.l. 
By breaking the three-fold degeneracy of the energy, this 
distortion selects a spin state with wavevector along x 
over its two twins with wavevectors rotated by ±7r/3. 

The amplitudes B21+1 and G21+1 are plotted versus 
A"i/|Ji| for D/lJxl = 0.1 in Fig.3(a). The phase <\> in 
Eq.([2]) is fixed by the phase of the lattice distortion. 
With increasing Ki/\Ji\, the third-order amplitude C3 
decreases while the amplitudes B\ and B 3 increase in 
size. For A'i/|Ji| = 0.05, B x w -0.21 and B 3 w -0.05. 
The anharmonicity of the spin changes with the distor- 



tion: (S z 2 ) decreases from 0.65 to 0.59 as K\ increases 
from to 0.05|Ji|. Spin population profiles P(S Z ) for 
A'i = and Ai = 0.05|Ji| are plotted in the inset to 
Fig. 3 (a). Their difference indicates that the distorted 
lattice is more heavily weighted towards S z = and the 
non-distorted lattice is more heavily weighted towards 
S z = ±1- As shown in Fig. 3(b), the |||| phase obtains 
a lower energy than the CNC phase when K\ > 0.057|Ji|. 

Other distortions of the lattice with a period of 1 
along x can also produce sizeable amplitudes B\ at the 
wavevector 27rx— Q. However, the 27rx — Q peak cannot 
be induced by the q = "isosceles" distortion observed 
by Feng et al. [22| in pure CuFe02 with the bonds 
along the x axis reduced in size but the diagonal bonds 
j} 2 ^ = J-[ 3 ' remaining identical. 

Guided by the observed elastic peaks at Q and 27rx— Q, 
Nakajima et al. Q constructed a CNC phase different 
than the one proposed here. Their "proper helix" is a 
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modified spiral with the same spin on sites R = mx + 
n-\/3y and R' = R + x/2 + V3y/2. After minimizing 
its energy as a function of wavevector for D/\Ji\ = 0.1, 
we compare the "proper helix" with the predicted CNC 
phase as well as with the pure 4444 phase in Fig. 3(b). Not 
only does the "proper helix" have a higher energy than 
the predicted CNC phase, it also has a higher energy than 
the |T 14 phase provided that D/\J\\ is not too small. For 
a non-distorted lattice with K\ = 0, the "proper helix" 
has a lower energy than the ||44 phase (but not lower 
than the predicted CNC phase) when D/\J X \ < 0.047. 

Estimating the actual distortion in CuFeC>2 requires 
that we also consider the elastic energy cost proportional 
to K\ 2 . Since the gain in exchange energy in Fig. 3(b) is 
linear in K\ , a distortion with period 1 occurs in both the 
4444 an d the CNC phases. Allowing a general distortion 
of the lattice with wavevector q, we find that the energy 
also has minima at q = 2Q and 47rx— 2Q, corresponding 
to a charge modulation with half the period of the spin 
modulation. This modulation has been observed Q in 
Al-doped CuFeC>2 and may be related to the predicted 
ferroelectric instability which contains both second- 
and fourth-order harmonics in addition to the uniform 
displacement of the oxygen atoms. 

Of course, it remains possible that the trial spin state 
used in this work is incomplete and that an even lower- 
energy state can be achieved. But the close agreement 
with Monte-Carlo simulations [H in Fig.3(b) bolsters 
our confidence that the anharmonic CNC state provides 
an excellent approximation to the true ground state of 
Eq.(p} for classical spins. As a test of our model, Fig. 3(a) 
indicates that the multifcrroic state should have small 
elastic peaks at the third harmonics of Q and Q' 

To summarize, we have shown that the CNC phase of 
a frustrated triangular lattice contains significant odd- 
order harmonics of the fundamental wavevector Q. This 
result should greatly facilitate the future modeling of 
multifcrroic ground states. As the easy-axis anisotropy 
D is reduced, the amplitudes of the higher harmonics are 
decreased and a pure cycloid is recovered as D — > 0. The 
27rx — Q peaks observed in Al- and Ga-doped CuFe02 
are explained by the scalene distortion of the triangular 
lattice. 
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